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Abstract 

We study K-type asteroids in the broad surroundings of the Eos family because they seem to be intimately related, according to 
their colours measured by the Sloan Digital Sky Survey. Such 'halos' of asteroid families have been rarely used as constraints for 
dynamical studies to date. We explain its origin as bodies escaping from the family 'core' due to the Yarkovsky semimajor-axis drift 
and interactions with gravitational resonances, mostly with the 9/4 mean-motion resonance with Jupiter at 3.03 AU. Our jV-body 
dynamical model allows us to independently est i mate t he age of the family 1 .5 to 1 .9 Gyr. This is approximately in agreement with 
the previous age estimate bv lVokrouhlickv et al.l (2006) based on a simplified model (which accounts only for changes of semimajor 
axis). We can also constrain the geometry of the disruption event which had to occur at the true anomaly / 120° to 180°. 
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1. Introduction 

The Eos family is one of the best-studied families in the main 
asteroid belt. Although we do not attempt t o repeat a thorough 
review presented in our previous paper IVokrouhlickv et al. 
(2006), we recall that the basic structure of the family is the 
following: (i) there is a sharp inner boundary coinciding with 
the 7/3 mean-motion resonance with Jupiter at approximately 
2.96 AU; (ii) the 9/4 mean-motion resonance with Jupiter di- 
vides the family at 3.03 AU and asteroids with larger sizes are 
less numerous at larger semimajor axes; (iii) there is an ex- 
tension of the family along the z\ = g - g(, + s - S(, secular 
resonance towards lower values of proper semimajor axis a p , 
eccentricity e p and inclination sin/ p . All these fact seem to be 
determined by the interaction between the orbits drifting due 
to the Yarkovsky effect in semimajor axis and the gravitational 
resonances which may affect eccentricities and inclinations. 

In this work, we focus on a 'halo' of asteroids around the 
nominal Eos family which is clearly visible in the Sloan Dig- 
ital Sky Survey, M oving Object Catalogue version 4 (SDSS, 
Parke r et al. 2008). As we shall see below, both the 'halo' and 
the family have the same SDSS colours and are thus most likely 
related to each other. Luckily, the Eos family seems to be spec- 
trally distinct in this part of the main bel t (several Eos fam- 
ily members were classified as K-types by DeM eo et al. 2009) 
and it falls in between S-complex and C/X-complex asteroids in 
terms of the SDSS colour indices. Detailed spectrosco pic ob- 
servations were also performed by IZappala et al.l (120001) which 
confirmed that asteroids are escaping from the Eos family due 
to the interaction with the J9/4 resonance. 

Our main motivation is to understand the origin of the whole 
halo and to explain its unusually large spread in eccentricity and 



inclination which is hard to reconcile with any reasonable ini- 
tial velo city field. Essentially, this is a substantial extension of 
work of Vokrouhli ckv et al. (2006), but here we are interested 
in bodies which escaped from the nominal family. 

We were also curious if such halos may be somehow related 
to the giant-planet migration which would have caused signif- 
icant gravitational pertu rbations of all small-body populations 
dMorbidelli et al. , 2005). Of course, in such a case the process 
is size-independent and moreover the age of the corresponding 
family would have to approach 3.9 Gyr in order to match the 
Nice model of giant-planet migration. 

In Section [2] we define the Eos halo and core populations. 
Section [3] is devoted to a description of our dynamical model 
and to a comparison with the SDSS observations. We discuss 
consequences of our results in Section [4] 

2. A discernment of the family core and halo 

In this Section, we proceed as follows: (i) we use a hierarchi- 
cal clustering method to extract the nominal Eos family; (ii) we 
look at the members of the family with SDSS colours and we 
define a colour range; (iii) we select all asteroids with Eos-like 
colours from the SDSS catalogue; finally, (iv) we define a halo 
and core using simple 'boxes' in the proper-element space. 

2.1. Colours of Eos-like asteroids 

We want to select asteroids similar to the Eos family, but 
first we have to choose a criterion to do so. We thus identify 
the no minal Eos family usi ng a hierarchical clustering method 
(HCM. IZappalaetalJ 1995) with a suitably low cut-off velocity 
^cutoff = 50 m/s (which leads to a similar extent of the family 



Preprint submitted to Icarus 



February 7, 2013 



0.4 



0.2 



-0.2 - 



-0.4 



1 .1 , _ 

* • -•*"■;"".** 

'y;.ijlajgjjgt 

V.''iaBBI 


1 . *■• 1 . 

- 1 • • • • ' - 




HHkjl' : ::. 


















1 1 1 





-0.4 



-0.2 0.2 
a I mag 



0.4 




-0.2 0.2 
a I mag 



Figure 1 : Colour indices i - z and a* (defined in Parker et al. ( 2008)) of all as- 
teroids from the Sloan Digital Sky Survey, Moving Object Catalogue version 4 
and the corresponding colour palette (top panel) which is used in the following 
figures to distinguish colours of asteroids. We also plot the Eos family mem- 
bers observed by the SDSS (bottom panel) with small photometric uncertain- 
ties (less than 0.03 mag). The inferred range of colour indices (denoted by the 
dashed yellow rectangle) is then used as a criterion for the selection of the Eos- 
like asteroids in the broad surroundings of the nominal family. The rectangle 
does not encompass the outliers. 



as in Vokrouhlickv et al and extract colour data from 

the SDSS catalogue (see Figure[T]i. The majority of Eos-family 
asteroids have colour indices in the following intervals 



a € (0.0, 0.1) mag, 
i-z e (-0.03, 0.08) mag 



(1) 

(2) 



which then serves as a criterion for the selection of Eos-like 
asteroids in the broad surroundings of the nominal family. 

We also used an independent method for the selection of Eos- 
like asteroi ds employing a 1 - dimensional colour index (which 
was used in Parker et aDl2008l to construct their colour palette) 



and we verified that our results are not sensitive to this proce- 
dure. 



Table 1 : The definitions of the core, halo and background populations in terms 
of intervals of proper eccentricity e p and proper inclination sin / p . The range of 
proper semimajor axis a p e (2.95, 3.16) AU is the same in all cases. 



population 


e 


p 


sin/p 


note 


core 


0.04- 


-0.10 


0.15-0.20 




halo 


0.00- 


-0.15 


0.12-0.24 


and not in the core 


background 


0.00- 


-0.15 


0.06-0.12 


together with. . . 




0.00- 


-0.15 


0.24-0.30 





while the range of proper semimajor axis is always the same, 
a p e (2.95,3.16) AU. 

Our results do not depend strongly on the selection criterion. 
For example, we tested a stringent definition: core was identi- 
fied by the HCM at v cut orF = 50 m/s and all remaining bodies in 
the surroundings belong to the halo. This approach makes the 
core as small as possible and the halo correspondingly larger 
but our results below (based on halo/core ratios) would be es- 
sentially the same. According to our tests, not even a different 
definition of the background/halo boundary changes our results. 

We are now ready to construct size-frequency distributions 
of individual populations. In order to convert absolute mag- 
nitudes H to diameters D we comput ed the median geom etric 
albedo p v = 0.16 from the WISE data (iMasiero et al .11201 lb for 
the nominal Eos family members. The size-frequency distribu- 
tion (Figure|4]i of the halo has a cumulative slope N(> D) cc D 7 
equal to y — -3.9 ± 0.2 in the size range D — 6 to 15 km and 
is significantly steeper than that of the core (y = -2.2 + 0.1). 
Even this difference of slopes ( 1 .7 + 0.2) indicates that if there a 
process transporting asteroids from the core to the halo it must 
be indeed size-dependent. 



A fr eq uency analysis s imilar as in lCarruba and Michtchenko 
(2007) or Carruba (2009) shows that there is approximately 5 % 
of likely z\ resonators (with the frequency g - g^ + s - S(, < 
0.3 " /yr) in the halo region. However, the concentration of ob- 
jects inside and outside the resonance is roughly the same, so 
that this secular resonance does not seem to be the most impor- 
tant transport mechanism. 



3. Yarkovsky-driven origin of the halo 

Motivated by the differences of the observed SFD's, we now 
want to test a hypothesis that the Eos family halo (or at least a 
part of it) was created by the Yarkovsky semimajor-axis drift, 
which pushes objects from the core into neighbouring mean- 
motion resonances and consequently to the halo region. 



2.2. Boundaries in the proper element space 

Next, we have to distinguish the family 'core' and 'halo' pop- 
ulations on the basis of proper orbital elements (a p ,e p ,sin/ p ) 
which will be consistently used for both the SDSS observa- 
tions and our dynamical models. We also need to define 'back- 
ground' population which enables to estimate how many as- 
teroids might have Eos-like colours by chance. We decided 
to use a simple box criterion (see Figures [2] [3] and Table [!}, 



3.1. Initial conditions 

We prepared an A^-body simulation of the long-term evolu- 
tion of the Eos core and halo with the following initial condi- 
tions: we included the Sun and the four giant planets on current 
orbits. We applied a standard barycentric correction to both 
massive objects and test particles to prevent a subst antial shift 



of secular frequencies (Mila ni and Knezevic , 1992). The to 



tal number of test particles was 6545, with sizes ranging from 
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Figure 2: The proper eccentricity e p vs proper inclination sin / p plot for as- 
teroids included in the SDSS MOC 4 catalogue. The proper semimajor axis is 
confined to the interval 2.95 to 3. 16 AU, i.e. the Eos family zone. Colour coding 
corresponds to the SDSS colour indices according to Figure The top panel 
includes all asteroids (regardless of their colours). The bottom panel shows 
only a subset of 'Eos-like' asteroids with colours similar to those of the Eos 
members (see Fig.[T] bottom). Moreover, we denote a box used for the defini- 
tion of the family 'core' (dashed yellow line) a larger box for the 'halo' (dotted 
green line) and two boxes considered as 'background' (thin black line). For 
comparison, we also plot positions of the nominal Eos family members (black 
dots), identified for the velocity v cuto ff = 50m/s. 
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Figure 3: The proper semimajor axis a p vs proper eccentricity e p (top panel) 
and Op vs proper inclination sin / p (bottom panel) for the observed Eos core 
(black dots), halo (red dots) and remaining Eos-like asteroids (gray dots) in the 
surroundings. The sizes of symbols are (inversely) proportional to the absolute 
magnitudes H of asteroids. The positions of important mean-motion resonances 
with Jupiter are also indicated (dotted vertical lines). 



10000 




10 

diameter Dl km 



100 



Figure 4: The cumulative size-frequency distributions N(>D) of the Eos core 
and halo. We show power-law fits and corresponding slopes y which clearly in- 
dicate that the halo population is significantly steeper than the core population. 
For comparison, we a lso plot the SFD of the nominal Eos family (as inferred 
from the WISE data, Masiero et al. 2011). The SFD's of the core and halo 
are biased because they include only asteroids observed by the SDSS. Conse- 
quently, the core seems to be much less populated than the nominal Eos family, 
even though these SFD's should be very similar. Nevertheless, the slopes and 
the halo/core ratios which we use in our analysis is not much affected by this 
bias. 
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D = 104 to 1 .5 km and the distribution resembling the observed 
SFD of the Eos family. 

Material properties were as follows: the bulk density p = 
2500 kg/m 3 , the surface density p s = 1500 kg/m 3 , the thermal 
conductivity K = 0.001 W/m/K, the specific thermal capacity 
C = 680J/kg/K, the Bond albedo A = 0.1, the infrared emis- 
sivity £ = 0.9, i.e. all typical values for regolith covered basaltic 
asteroids. 

Initial rotation periods were distributed uniformly on the in- 
terval 2 to 10 hours and we used random (isotropic) orientations 
of the spin axes. The YOR P model of t he spin evolution was 
described in detail in iBroz et al.1 (1201 II) . while the efficiency 
parameter was cy orp = 0.33 (i.e. a likely value according to 
Hanus et alj|201 ll) . YORP angular m omenta affecting the spin 
rate a nd the obliquity were taken from lCapek and Vokrouhlickv 
(2004). We also included spin axis reori entations caused by col - 
lisions^ with a time scale estimated by Farinella et al. (1998): 
Tteot = Bitu/cjof 1 (D/Dof 2 , where B = 84.5 kyr, B x = 5/6, 
Bi = 4/3, Do = 2 m and a»o corresponds to period P — 5 hours. 

The initial velocity field was size-dependent, v oc vqDq/D, 
with vo = 93m/s and Dp = 5 km (i.e. the best-fit val- 
ues from IVokrouhlickv et al.ll2.Q06h . In principle , this type of 



size-v elocity relation was initially suggested by Celli ng et al 



(1999), but here, we attempt to interpret the structure of the 
family as a complex interplay between the velocity field and 
the Yarkovsky drift which is also inversely proportional to size. 
We assumed isotropic orientations of the velocity vectors. The 
geometry of collisional disruption was determined by the true 
anomaly / = 150°, and the argument of perihelion cj - 30°. 
We discuss different geometries in Section |4] 

We use a modif ied version of the SWIFT package 

( Levison and Duncan! 1994 ) for n umerical integrations, with 
a second-order symplectic scheme dLaskar and Robutefl 2001 ), 
digital filters employing frequency -modified Fourier transform 



( Sidl ichovskv and N esvornv, 1996) and an implementation of 



the Yarkovsky effect dBrozl 12006b . The integration time step 
was At = 91 days, the output time step after all filtering proce- 
dures 10 Myr and the total integration time span reached 4 Gyr. 

3.2. Results of the N-body simulation 

Initially, almost all asteroids are located in the core (see Fig- 
ure |5). Only a few outliers may have velocities large enough to 
belong to the halo. Within a few million years the halo/core 
ratio quickly increases due to objects located inside the 9/4 
resonance and injected to the halo by these size-independent 
gravitational perturbations. Further increase is caused by the 
Yarkovsky/YORP semimajor axis drift which pushes additional 
orbits into the J9/4 and also other resonances. 

We checked the orbital elements of bodies at the moment 
when they enter the halo region (Figure [6} and we computed 
the statistics of dynamical routes that had injected bodies in the 



We do not take into account collisional disruptions because we model only 
that subset of asteroids which survived subsequent collisional grinding (and 
compare it to the currently observed asteroids). Of course, if we would like 
to discuss e.g. the size of the parent body, it would be necessary to model 
disruptive collisions too. 
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Figure 6: The proper semimajor axis a p vs proper inclination sin/ p of the syn- 
thetic family members at the moment when they have entered the halo region. It 
is easy to distinguish objects injected by mean-motion resonances and by sec- 
ular resonances in this projection, because the former have a particular value 
of the semimajor axis. The objects injected by the J9/4 resonance are denoted 
by red colour, the Jl 1/5 and 3J - 2S - 1 by orange, the J7/3 by yellow, the Z[ 
secular resonance and other resonances by black. 



halo: J9/4 57%, Jl 1/5 (together with a three-body resonance 
3J - 2S - 1 with Jupiter and Saturn) 10 %, J7/3 6 %, and z\ sec- 
ular resonance 23 %. The remaining few percent of bodies may 
enter the halo by different dynamical routes^ However, if we 
account for the fact that bodies captured by the z\ resonance 
usually encounter also the J9/4 resonance that scatters them fur- 
ther away in to the halo, we obtain a modified statistics: J9/4 
70 %, Jl 1/5 12 %, J7/3 5 %, and zi 10 % that better reflects the 
importance of different mechanisms. 

A saturation of the halo occurs after approximately 1 Gyr, 
because the halo population is affected by the Yarkovsky/YORP 
drift too, so that the injection rate roughly matches the removal 
rate. Nevertheless, the halo/core ratio steadily increases, which 
is caused by the ongoing decay of the core population. 

In order to compare our model and the SDSS observations 
we compute the ratio % = dAfhaio/dAf core between the number 
of objects in the halo and in the core for a given differential 
size bin. This can be computed straightforwardly from our sim- 
ulation data. In case of the SDSS observations, however, we 
think that there is a real background of asteroids with Eos-like 
colours (may be due to observational uncertainties or a natural 
spread of colours; see Figure [2}. Obviously, such background 
overlaps with the core and the halo, so we need to subtract this 
contamination 



^obs = 



dA^halo - 0.833 (Wbackgroun 
dJVcore - 0.167 dA^backeioun 



(3) 



The numerical coefficients then reflect different 'volumes' of 
the halo, core and background in the space of proper elements 
(a p , e p , sin/p), as defined in TableQ] 

As we can see in Figure|7] a reasonable match to the observed 
halo/core ratios can be obtained for ages 1.5 Gyr (for smaller 



2 Other secular resonances intersecting this region, s — s$ - 2g$ + 2g^ or 
g + 2g$ - 3g(,, do not seem to be important with respect to the transport from 
the core to the halo. 
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Figure 5: Left panels: the proper semimajor axis a p vs proper eccentricity e p plots showing a dynamical evolution of our synthetic family. We can distinguish the 
core (black dots), the halo (red dots) and objects beyond the halo box (gray dots). There is a comparison to the observed Eos core and halo too (yellow crosses), 
as inferred from the SDSS data (the same as in Figure [2] bottom). The positions of important resonances are indicated by vertical dotted lines. We plot the initial 
situation at t = (top panel) and the evolved family at t = 1.7 Gyr (bottom panel). The core of the synthetic family exhibits a slightly different structure than the 
observed core which may indicate that: (i) the initial true anomaly was closer to / = 180°, or (ii) the initial velocity field deviated from the assumed v ct I/O 
dependence. Right panels: the corresponding size-frequency distributions of the synthetic core (black line), which was always scaled to the observed SFD of the 
Eos core, and the synthetic halo (red line) which can be then directly compared to the observed halo (gray line). It is clear that the halo's SFD becomes steeper in 
the course of time and at t = 1.7 Gyr it matches the observed SFD. 
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Figure 7: Top panel: the evolution of the halo/core ratio in our simulation for 
several size bins (colour curves) and its comparison to the observed SDSS ra- 
tios in the same bins (horizontal lines). The intersections give age estimates 
from 1.5 Gyr (for smaller bodies) to 2.2 Gyr (larger bodies). Bottom panel: the 
corresponding evolution of the x 1 vs ti me t- The dotted line indicates the num- 
ber « = 8 of size bins (from D = 2 km to 10 km) in which the^ 2 was computed. 
The best fits (with^ - ss n or smaller) correspond to ages from 1.5 to 2.2 Gyr. 
For comparison, we also plot a^ 2 (r) dependence (red line), computed from the 
histog ram of C = (a - 3.019 AlT l/lO^ values, which corresponds to the anal- 
ysis of Vokrouhlickv et al. ( 2006). We use the population of Eos-like asteroids 
observed by the SDSS in both core and halo for this purpose. 



bodies) to 2.2 Gyr (for larger bodies). To better quantify the 
difference between the model and the observations we construct 
a suitable metric 



X 



i=2 



(X(0-Kbs,) 2 



(4) 



where the summation is over the respective size bins (£),-, D, + 
dD), Dj = i ■ 1 km and AD = 1 km. The uncertainties of the 
numbers of objects are of the order <x halo VdA^aio, 0~core - 
VdA^ core , and cr, reflects their propagation during the calculation 
of the ratio in Eq. Q in a standard way 



O"/ = -yVhalo/dNhalo) 2 + (0"core/dA^ core ) 2 djVhalo/dN core 

and similarly for <x b SI . The x 2 (t) dependence is shown in 
Figure [7] and the best-fit is obtained again for the ages t 
1.5 to 2.2 Gyr. 

The ratios % are directly related to the size-frequency dis- 
tributions and consequently we are indeed able to match the 
observed SFD's of halo and core, including their slopes and ab- 
solute numbers (Figure [5] right column) 

These results are not very sensitive to the initial velocity 
field, because most asteroids fall within the family core; ve- 
locities would be unreasonably large to have a substantial halo 
population initially. 



4. Conclusions 

Yarkovsky-driven origin seems to be a natural explanation of 
the halo population. A lucky coincidence that the disruption 
of the Eos-family parent body occurred close to the moderately 
strong 9/4 mean motion resonance with Jupiter established a 
mechanism, in which orbits drifting in semimajor axis due to 
the Yarkovsky effect are mostly perturbed by this resonance 
and scattered around in eccentricity and inclination. The total 
spread of the simulated halo (up to 0.2 in eccentricity, Figure[5), 
which matches the SDSS observations (Figure|2), also supports 
our conclusion. 

As an important by-product, the process enabled us to in- 
dependently constrain the age of the family. Moreover, if we 
analyse the evolution in the proper semimajor axis vs the ab- 
solute magnitude (a p , H) plane and create a histogram of the 
quantity C = (a - 3 . 019 A U)/10 W (i.e. a similar approach as in 
IVokr ouhlickv et al.1 d2006h . but now using a full A^-body model 
and the SDSS observations for both the core and halo), we can 
compute an independent ^ 2 (f) evolution (refer to Figure Q red 
line). Since both methods - the halo/core ratios % and the C- 
histogram - seem to be reasonable, we can infer the most prob- 
able age as an overlap of intervals of low x 2 (t) and this way 
further decrease its uncertainty, so that t 1 .5 to 1 .9 Gyr. 

It is also interesting that the true anomaly at the time of dis- 
ruption has to be / a 120° to 180°. We performed tests with 
lower values of / and in these cases the synthetic family has 
initially a different orientation in the (a p , e p ) plane: the objects 
are spread from small a p and <? p to large a p and e p (cf. Figure|5]l. 
Way too many objects thus initially fall in to the z\ secular reso- 
nance and because such captured orbits cannot drift to small a p 
and large e p it is then impossible to explain the observed struc- 
ture of the family and consequently / < 120° is excluded. 

Finally, let us emphasize that given the differences between 
the size-frequency distribution of the halo that of the core, we 
can exclude a possibility that the Eos halo was created by a 
purely gravitational process (like the perturbations arising from 
giant-planet migration). 
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